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Abstract: A modification of the relation between axion mass and the PQ con- 
stant permits a relaxation of the astrophysical constraints, considerably enlarging 
the allowed axion parameter space. We develop this idea in this paper, discussing 
a model for an ultramassive axion, which essentially represents a supersymmetric 
Weinberg- Wilczek axion of the mirror world. The experimental and astrophysical 
limits allow a PQ scale f a ~ 10 4 — 10 6 GeV and a mass m a ~ MeV, which can be ac- 
cessible for future experiments. On a phenomenological ground, such an ultramassive 
axion turns out to be quite interesting. It can be produced during the gravitational 
collapse or during the merging of two compact objects, and its subsequent decay into 
e + e~ provides an efficient mechanism for the transfer of the gravitational energy of 
the collapsing system to the electron-positron plasma. This could resolve the energy 
budget problem in the Gamma Ray Bursts and also help in understanding the SN 
type II explosion phenomena. 
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1. Introduction 

The strong CP problem is one of the most puzzling points of modern particle physics 
(for a general reference see e.g. Jl|). It resides in the presence of the so-called 9— term 
in the QCD Lagrangian, C e = 9(a s /8n)G ^ U G^ V ', where a s is the fine structure con- 
stant of the strong interactions and is the gluon field strength tensor. The 
9 -term, which receives a contribution also from the complex phases in the quark 
mass matrices M UjD so that its effective value is 9 = 9 + arg(det Mu det Md), is CP 
violating and leads to a neutron electric dipole moment d n , experimentally not ob- 
served. This implies a very strong upper limit on the parameter 9, \9\ < 10~ 9 , which 
has no theoretical explanation in the context of the Standard Model. 

In the most appealing solution of this problem, the Peccei-Quinn (PQ) mech- 
anism ||, the 9 parameter 1 becomes a dynamical field, the axion a = f a 9, and 
emerges as the pseudo-Goldstone mode of a spontaneously broken global axial sym- 
metry C/(1)pq. Here /„ is a constant, with dimension of energy, also called axion 
decay constant. We will use the following convention throughout this paper |J: we 
indicate with /pq the VEV of a scalar (or a VEV of a combination of several scalars) 
responsible for the C/(1)pq symmetry breaking. The constant /„, which characterizes 
axion phenomenology, is defined as /pq/iV, where iV stands for the color anomaly of 
£7(1) pq current. 2 

1 In the following we will use 9 instead of 9 for simplicity. 

2 The PQ charges are normalized so that each of the standard fermion families contributes as 
N = 1. Therefore, in the Weinberg- Wilczek (WW) model || we have N = N g , where N g (= 3) is 
the number of fermion families. The same holds in the Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) 
model Q. Other models of the invisible axion, e.g. the hadronic axion (f| or archion ||, generally 
contain some exotic fermions and so N ^ N g . 
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At quantum level, the PQ symmetry is broken by the chiral anomaly, and since 
this is a dynamical effect, its strength is measured by non-perturbative QCD con- 
tributions. Therefore all the axion properties are essentially related only to the PQ 
scale /pq and the QCD scale A. 

On a phenomenological ground, all the axion characteristics can be roughly es- 
timated from the pion properties. Indeed, axions generically mix with pions so that 
their mass, as well as their couplings with photons and nucleons, are roughly given 
by /„•/ f a times those of the 7r-meson, where f n ~ 93 MeV is the pion decay constant. 
So it is clear that the PQ constant is the relevant scale for axion phenomenology. For 
example, in the most general context, axion interaction with fermions and photons 
can be described by the following Lagrangian terms: 



C-a = iCae-T fl ^ G + ic aN^~ A75A + C ay — — a F^F^ H , (1.1) 

Ja la OTTfa 



where represents the fermion mass (e.g. m e ,mN, ... for electrons, nucleons, etc.), 
a is the fine structure constant and c ai are constant coefficients. The factors c a N, 
which refer to axion-nucleon interaction, N = p,n, are generally of order one, while 
for the axion-electron and axion-photon interaction the related coefficients c ae , c a7 
are model dependent. 

As for the axion mass, in general it can be obtained from the expression : 

2 _ 1 VK 
ma ~ tfV + KTrM-i ' (L2) 

where V = (qq) ~ A 3 is the light quark condensate responsible for the chiral sym- 
metry breaking, M = diag(m u , rrid, •••) is the mass matrix of light quarks (with 
m q < A) and K ~ A 4 accounts for the strength of the instanton induced potential. 
Then, by taking into acount only two light quarks, u and d, and using the relation 
(m u + m ( j)(qq) = Ta^f^, from ( |1.2j ) one directly arrives to the more familiar formula: 

m = ]_ f mumgV X 1 ' 2 = _ K ^ AO 6 GeV \ x g 2 gV (1 3) 

fa\m u + m d J l+Z f a \ f a J ' ' 

where z = m u /m a ' ~ 0.57. 

In the WW model the PQ symmetry is broken by two Higgs doublets 
with the VEVs v 1>2 , namely f a = (v/2) sin 2/5, where v — (vf + f|) 1/2 ~ 247 GeV is 
the electroweak scale and tan/3 = Vijv\. Therefore, the WW axion is quite heavy, 
and its mass can vary from a few hundred keV to several MeV: 

ww 2 A fm u V\ 1/2 _ 150 keV 



WW = ~ _ _ ■ (l A) 

vsin2/? \l + z) sin2/5 ' 1 ' ; 

However, its couplings with fermions and photons are too strong and for this reason 
the WW model is completely ruled out for any values of the parameter (3 by a variety 
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of terrestrial experiments such as the search of the decay K + — > n + a, the J/ip and 
T decays into a + 7, the nuclear deexcitations via axion emission, the reactor and 
beam dump experiments, etc. ||. 

For any realistic axion model, these experimental constraints generally imply the 
lower bound on the PQ scale, f a > 10 4 GeV, which in turn implies m a ^ 1 keV. This 
is what happens in all the so-called invisible axion models. Anyway, such a light 
axion can easily be produced inside a star at a temperature of a few keV, and can 
accelerate the cooling process in a dangerous way. So the axion is required to interact 
only weakly with fermions in order to strongly suppress the energy transport process 
inside the star. Also, the axion luminosity from the SN core must be constrained in 
order to not ruin the neutrino signal detected at the time of SN 1987A explosion [|]. 
These astrophysical considerations exclude all scales f a up to 10 10 GeV [[I], and so 
m a < 10 -4 eV. 3 On the other hand, the cosmological limits related to the primordial 
oscillations of the axion field or to the non-thermal axion production by cosmic 
strings, demand the upper bound f a < 10 11 — 10 12 GeV [[l], ||. Thus, not much 
parameter space remains available. 

Is it possible to relax the astrophysical constraints? It is quite interesting to note 
that all astrophysical constraints from stellar evolution could be satisfied for the PQ 
scales above the laboratory limit, f a > 10 4 GeV, up to values order 10 7 GeV, if the 
axion mass would remain in the MeV range, as with the mass of the WW axion ( |1.4j ). 
However, the tight relation between the axion mass m a and the PQ scale f a ( |1.3|) is 
very constraining and does not allow such a situation. 

If there were, in fact, another source for axion mass, e.g. from Plank scale induced 
effects, then the axion could change its properties and then no longer be useful for 
the solution of the strong CP problem. This explains why this relation is universally 
accepted, and in several papers the astrophysical and cosmological limits are given 
in terms of axion mass instead of the PQ constant. 

But is the relation ( |1.3|) really universal? What would happen, in fact, if the 
axion could communicate with another, hidden, sector of particles and interactions? 
In general, if the axion potential were to get dominant contribution from the hidden 
sector, then we would expect it to solve the strong CP problem for that sector rather 
than for our observable world. However, this is not necessarily the case. 

In particular, one can assume that the hidden sector is a mirror world, a parallel 
sector of "mirror" particles and interactions with the Lagrangian completely similar 
to that of the ordinary particles [fLOfl . In other words, it has the same gauge group 
and coupling constants as the ordinary sector, so that the Lagrangian of the whole 
theory is invariant with respect to the Mirror parity (M-parity) under the interchange 
of the two sectors. Several phenomenological and astrophysical implications of the 



3 In the case of the hadronic axion [|| or archion |6j, a small window around f a ~ 10 6 GeV (axion 
mass of a few eV) can be also permitted. 
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mirror world have been studied in the literature [11|. In particlular, it could provide 
a new insight for the cosmological dark matter made up by mirror baryons |T^| . For 
a recent review of the mirror matter concept one can refer to [U 



We can further assume that the ordinary and mirror sectors have the common 



Peccei-Quinn symmetry [1^], with the same PQ charges carried by the ordinary Higgs 
doublets Hi 2 and their mirror partners H[ 2 . 

If the M-parity is an exact symmetry, then the particle physics should be exactly 
the same in two worlds, and so the strong CP problems would be simultaneously 
solved in both sectors. In particular, the axion would couple to both sectors in 
the same way and their non-perturbative QCD dynamics would produce the same 
contribution to the axion effective potential. This situation would not bring drastic 
changes of the axion properties; just the axion mass would increase by factor of \/2 
with respect to the standard expression (|1.3|). 4 

However, a very interesting situation emerges if the M-parity is spontaneously 
broken in the Higgs sector, so that the mirror electroweak scale v' is considerably 



larger than the ordinary one v |TB|, |T7]. In this case, as it was shown in ref. (T7|, |T4| ], 
also the mirror QCD scale A' becomes larger than the ordinary one A and thus one 
expects that the dominant contribution to the axion potential comes from the mirror 
sector. On the other hand, if the latter contribution is predominant, the axion mass 
relation with the PQ scale could change in a considerable way. Of course, in this 
case it is absolutely not evident that the axion can still solve the strong CP problem. 
Indeed, if we ask for it to be dominantly governed by the mirror QCD, it is natural 
to expect that it will cancel the mirror and not the ordinary ^-term. Fortunately, 
this is not the case: this axion can still solve our strong CP problem, as far as we 
ask for the Yukawa structure to be the same in the two sectors [fL4| . So it appears 
that the mirror extension of the standard axion model is the only way that allows 
for a quite heavy axion without spoiling the PQ mechanism. 

In particular in this paper we will consider an axion with the mass m a > 1 MeV 
and with a PQ scale f a ~ 10 6 GeV. Certainly such a massive axion cannot ruin 
stellar evolution process, even though it can still be produced in the hot SN core at 
a temperature of a few 10 MeV. We expect then that the available parameter range 
is sizeably increased with respect to the standard DFSZ axion model. 

As we will show, such an axion can have many interesting astrophysical impli- 
cations. In particular, it can help in understanding the Gamma Ray Bursts (GRBs) 
and the supernova explosion phenomena via the mechanism suggested in ref. |19| 6 . 



4 However, even in this less interesting case there would be significant modification of the pri- 
mordial oscillations of axion field and their contribution to dark matter of the Universe |L5| . 

5 Cosmological implications of such a scenario were studied in refs. [jl7|, |l8). 

6 We underline that the possible relation between GRBs and supernova was contemplated in the 
past. In pOj , e.g., is considered the possibility that the emission of a light axion from SN can 
provide the energy necessary for the GRBs production. 
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In fact, such a non standard axion can be produced during the merging of two com- 
pact objects, and then, due to its decay, efficiently transfer the gravitational energy 
of the collapsing system into an ultrarelativistic e + e~ plasma, the fireball, far from 
the impact place. Also such a heavy axion, produced in the core of a SN type II, 
can decay into e + e~ before reaching the stellar surface, transferring in this way a 
huge amount of energy at a distance of about 1000 Km from the stellar core, helping 
the SN explosion (thermal bomb). On the other hand, because of the different size 
of SN type Ib/c, some axions will be able to leave their surface and then decay into 
photons, explaining the observed events of weak GRBs related to the SN type I. 

This logic was essentially suggested in a previous paper [Q, but the allowed 
mass was not enough to open a new window for the axion of interest in astrophysics. 
In fact, the parameter space in |14[ does not allow for an f a > few 10 4 GeV. Observe 
that the QCD scale in the mirror sector (and , as a consequence, the axion mass) is 
fixed by the renormalization group (RG) evolution of the strong coupling constant, 
and so it depends on the matter content of the theory, as well as on the mirror 
fermion masses. This observation suggests to us that a supersymmetric version of 



the model can lead to a quite different phenomenology. The results obtained in |14 



are summarized in figure 2 of the cited paper. We observe that the described axion 
is constrained, from SN data, in the region between the dot-dashed line and the 
dashed one. Also, from terrestrial experiments, the PQ constant f a must be bigger 
than 10 4 GeV if the axion mass stays below ~ 1 MeV (twice the electron mass), 
otherwise the limit on f a should be increased to f a > few 10 5 GeV. We see that 
for m a > 1 MeV, the raising of the dashed line is strongly suppressed. This is 
because, for an axion heavier than two electrons, the total axion decay width T tot is 
dominated by the axion decay into electrons and positrons. Thus the photon flux is 
strongly suppressed. This phenomena can be improved if the axion mass is increased 
for a given PQ constant f a , since T(a — > e + e~) ~ m a / f%. This, as we will show, 
is exactly what happens in a supersymmetric model, for a certain range of f a . We 
will then consider a supersymmetric mirror axion model and this will allow us to 
strongly enlarge the parameter space for the axion with relevant phenomenological 
consequences. 

Finally, in this model the hierarchy problem between the mirror and the ordinary 
e.w. scale is solved naturally via the GIFT (Goldstones Instead of Fine Tuning) 
mechanism [BT|. In fact, the Higgs potential has an additional SU{4) symmetry. This 
accidental symmetry is global, but contains the local one SU{2) x SU{2) = 5*0(4). 
When the M-parity is broken, the mirror Higgses acquire a large VEV, while the 
ordinary ones remain as Goldstone bosons until the supersymmetry breaking is taken 
into account. The ordinary Higgs VEVs are then naturally of the order of the 



supersymmetry soft breaking scale [21, 22 



The paper is organized as follows. In sect. 2, we present a simple mirror super- 
symmetric axion model and describe the general features and the experimental and 
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astrophysical bounds, showing that there is an allowed parameter space which can 
be of interest for the future experimental search. In sect. 3, we carefully study the 
interesting relations with the physics of the GRBs and with the dynamics of different 
supernovae. Finally in sect. 4, we summarize our results. 



2. The Mirror Weinberg- Wilczek axion 

Let us consider a model based on the gauge symmetry G x G' where G = SU (3) x 
577(2) x U(l) stands for the standard model of the ordinary particles: the quark 
and lepton fields ipi = q iy k, u\, d\, e\ [i is a family index) and two Higgs doublets 
Hi j2 , while G' = SU(3)' x SU(2)' x U(l)' is its mirror gauge counterpart with the 
analogous particle content: the fermions ^ = q[, l'^ uf, df, ef and the Higgses 



H[ 2 |10[ . From now on, all fields and quantities of the mirror sector have an apex to 
distinguish them from the ones belonging to the ordinary world. All fermion fields 
are taken in a left-chiral basis. 
Let us assume that the theory is invariant under the mirror parity M: G <-> G", 
which interchanges all corresponding representations of G and G' . Therefore, the 
two sectors are described by identical Lagrangians and all coupling constants (gauge, 
Yukawa, Higgs) have the same pattern in both of them. In particular, for the Yukawa 
couplings 

£ Yuk = Gg'ufo# 2 + G%d\ qj H x + G^etljHi + h.c, 

£ Yuk = G'ffufq'^ + G'iJdW 3 H[ + G'iiefl'H', + h.c, (2.1) 

we have G)j D E = G'ff D E . In addition, the initial ^-terms are equal, 6 = 9'. 

We further assume that the two sectors have a common Peccei-Quinn symmetry 
£/(1)pq realized a la Weiberg- Wilczek model. The essential point is then to have a 
term in the Higgs potential which mixes the Higgses of different kinds. The simplest 
possibility is: 

V mix = -KfaHrfiHiHti + h.c, (2.2) 

where the coupling constant k should be real due to M-parity. 7 

As far as the M-parity is an exact symmetry, the particle physics should be 
exactly the same in the two worlds: in particular, the quark mass matrices are 
identical, Mjj,d = M' UD) the QCD scales coincide, A = A', and the axion couples 
with both QCD sectors in the same way: f ( ^ 1 a(GG + G'G'), so that their non- 
perturbative dynamics should produce the same contributions to the axion effective 
potential. Clearly, in this case, the strong CP problem is simultaneously solved in 



7 Notice that in the limit k = there emerge two separate axial global symmetries, U{1)a for 
the ordinary sector under which ipi — ► exp(— iu>/2)tpi and — » exp(iuj)H 1,2, and U(\)' A for the 
mirror sector: ip^ — > exp(— ioj' /2)ip' i and H [ 2 — > exp(ia/)7J{ 2 . Therefore, the term V m i x demands 
that u>' = cu and thus it reduces U(1)a x U(1)' a to its diagonal subgroup t/(l)pQ. 
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Figure 1: Example of the renormalization group evolution of the strong coupling constants 
a s and a' s (respectively solid and dashed), as a function of the energy scale /x, for v' jv = 
10000. Super symmetry is supposed to be broken at the scale m s = mi. 



both worlds - the axion VEV cancels the #-terms both in the mirror and ordinary 
sectors. In such a realization, however, f a remains order 100 GeV and thus it is 
excluded on the same phenomenological grounds as the original WW model. 

The situation is more interesting when the M-parity is spontaneously broken, as 
it was suggested in ref. [[16], and the electroweak symmetry breaking scale v' in the 
mirror sector becomes larger than the ordinary one v. Since [7(1)pq is a common 
symmetry for the two sectors, the PQ scale is determined by the larger VEV v', i.e. 
fpQ ~ (i//2) sin2/3'. The axion state a dominantly comes from the mirror Higgs 
doublets H[ 2 , up to small (~ v/v') admixtures from the ordinary Higgses Hi 2 - 
Hence, it is a WW-like axion with respect to the mirror sector, while it couples with 
the ordinary matter as an invisible DFSZ-like axion. This would lead to somewhat 
different particle physics in the mirror sector and it is not a priori clear that the 
strong CP problem can still be simultaneously fixed in both sectors. However, as 
it was shown in [TJ]], this is just the long as the Yukawa structure in the 

two sectors is the same. It happens then that the mirror quark masses are scaled 
linearly with respect to the ordinary ones: m' uct = (2m u ,c,t, m 'dsb = d m d,s,b, where 
C2 = v' 2 /v2 = C(sin/3 / / sin/3) and £1 = v[/vi = ((cos j3'/ cos/3). At very high energies, 
fi ^> v', the strong coupling constants a s (/i) and a' s (jj) should be equal due to M- 
parity. Under the renormalization group (RG) evolution they both evolve down in 
parallel ways until the energy reaches the value of mirror-top mass m' t ~ £2^- Below 
it, a' s will have a different slope than a s , and this slope will change every time below 
the mirror quark thresholds m' b ~ Ci m b, etc. In the evolution of a s these thresholds 
occur at lower scales, fi = m t , mj, etc. Then it is very easy to determine the scale A' 
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at which a' s becomes large, once we know that for the ordinary QCD this happens 
at A ~ 200 MeV. In other words, A' becomes a function of £1,2, and for v' 3> v one 
could obtain a significant difference between the QCD scales: A' > A (see figure [I]). 

The value of the QCD constant in the mirror sector depends on the RG evolution 
of the mirror strong coupling. This, in turn, depends on the matter content of the 
theory. Therefore in a supersymmetric theory we expect a different result. In general, 
the relation between the ordinary and mirror QCD scale can be written as 

^ = A(ft,ft')(», (2.3) 

where A is a function of the angles ft and ft', while p is a constant. They both depend 
on the number of mirror light quarks. We can estimate, 8 for a non supersymmetric 
theory, p ~ 0.36, while in a supersymmetric theory we find the bigger value 9 p ~ 0.54. 

In the following we will consider the renormalizable Lagrangian described by the 
superpotential: 

W = T(yn - M 2 ) + S(HiH 2 + H[h' 2 - 7l 2 ), (2.4) 

where M is an energy scale M 3> v, Hi,H2,H'i,H 2 are the Higgs superfields and 
T , y, 1Z, S, are other superfields. 71, y, S have respectively the PQ charges uj-ji = 
—ujy = 1/2, us = —1, while for the Higgs fields uii = sm 2 ft, u 2 = cos 2 ft and 
the same for u[ and u>' 2 with ft replaced by ft'. T does not transform under PQ 
symmetry. For convenience we indicate with T, Y, R, S, Hi, H[ the scalar components 
of the superfields T, y, 1Z, S, Hi, H[. The second term on the right hand side of ( |2.4| ) 
fixes the VEV pattern of the ordinary and mirror Higgses. The constraint equation 
is (HiH 2 ) + (H[H 2 ) ~ (R 2 ) where (R 2 ) is fixed in the first term on the right hand 
side of ( |2.4j ): (R) ~ (Y) ~ M. Since the ordinary electroweak scale (H) is of the 
same order as the supersymmetry breaking scale m s , the previous equation fixes 
(H') ~ (R) ~ M. Therefore the R VEV fixes the M-parity breaking scale of the 
theory. Since M ^> v ~ m s we expect all the mirror fermion masses to be heavier 
than the ordinary ones and, as a consequence, A' > A. 

The complete M-invariant scalar potential consists in the sum of the F and D 
terms, plus the supersymmetry breaking contribution: 

V = V F + V D + V B . (2.5) 

More explicitly 

V F = V(H, H') + \TR\ 2 + \YR - M 2 \ 2 + \TY - 2RS\ 2 , (2.6) 

8 We assume, for simplicity, that there are no light quarks in the mirror sector. This is easily 
verified when f a ^ few 10 4 GeV. 

9 For this estimation we have assumed that the supersymmetry is broken at the top quark mass. 
It is also assumed that there are no mirror light quarks, and that the lightest mirror quark is heavier 
than the ordinary top quark. These assumptions are verified for f a > few 10 6 GeV. 
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where 



V(H, H') = \H X H 2 + H[H' 2 -R 2 \ 2 + S 2 {Hl + H 2 2 + H' 2 + H' 2 ) + h.c, (2.7) 
while: 



v - g2 



{HiH 2 ) 2 + (H' 1 + H 2 ) 2 + {2.1 



where (7 is the coupling constant for isospin and g y refers to the hypercharge. Finally 
Vb contains all the possible soft supersymmetry breaking terms and fixes the scale 
of the supersymmetry breaking m s . 

In this model the hierarchy problem between the mirror and ordinary Higgs VEVs 
is solved in a rather natural way by the Pseudo-Goldstone or GIFT mechanism |H] . 
Namely, while the mirror Higgses get VEVs order M 3> Mw, the ordinary Higgses 
can get the masses (and hence the VEVs) of the order m s . 

More precisely, in the supersymmetric limit Vb = 0, the F-term potential ( |2.7| ) 
has an accidental global symmetry SU(4), larger than the local symmetry SU(2) x 
SU(2)' acting on the Higgses H\ 2 and H[ 2 respectively. Therefore, if the non-zero 
VEVs are located on the mirror Higgses, then we obtain v[v' 2 = R 2 ~ M 2 , and the 
D-term Q2.8p gives tan/?' = v' 2 /v[ = 1, while the ordinary Higgses H\ t 2 remain as 
Goldstone superfields. 

Then, considering the soft terms V#, one generates the soft mass terms for H 



1.2 



and, as a remarkable property of the GIFT mechanism [21] , |22 ] , also supersymmertic 



/x-term, \iR\R 2 with fi ~ m s , as far the VEV (S) ~ m s is generated by the soft 
terms Vb- As a result, one can generate non zero VEVs V\,v 2 ~ m s (with v \ + v 2 = 
v 2 = (247 GeV) 2 ), and the parameter tan/? = v 2 /vi in general can be different from 
1. On the other hand, given that the mirror VEVs are generated at a large scale, 
v i = v 2 ~ the soft terms Vb cannot significantly shift their values and thus we 
remain with tan/3' w I. 

Let us discuss now in more detail the axion phenomenology in this model. The 



PQ scale is (see (I|) /pq = yj^jivj^j) 2 where Uj and Vj are respectively the PQ 
charges and VEVs of the scalar fields. We find / PQ = (f 2 + f 2 + (S) 2 + 1/2M 2 ) 1 / 2 
where / = (f/2)sin2/5 and /' = (v'/2) sin2/3'. In the interesting physical limit 
/' ~ M ^> / ~ m s we see that /pq is essentially equal to /'. 

For what concerns the axion mass, its square is in general given by the sum 
of the right hand side of ( |1.2|) and a similar term, but with V, K, M replaced by 
V, K', M' . The meaning of these parameters is obvious, M'is the mass matrix of 
the mirror light quarks, and the values K' ~ A' 4 and V ~ A' 3 characterize the mirror 
gluon and quark condensates. Assuming, for simplicity, that there are no light quarks 
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Figure 2: Pictorial representation of the mirror QCD scale and mirror axion mass, in 
relation to the ordinary value. All the three lines correspond to tan 0' = 1 and, from the 
bottom, to tan/3 = 5, tan/3 = 25, tan/3 = 55. Supersymmetry is supposed to be broken at 
the scale m s = rat. 



in the mirror sector and that A' ^> A, we find m a ~ -4—. The axion mass is then 

J a 

essentially driven by the mirror QCD. This mass can be several times bigger than 
the ordinary WW mass: 

ma -***PI A '\ C = 2(, m \ A ) 1/2 ~0.2. (2.9) 




m ww C \A J ' \{l + z)A / 



For the numerical computation we have used the following values of the quark masses: 
m u = 4 MeV, m d = 7 MeV, m s = 150 MeV (at fj, = 1 GeV), and m c (m c ) = 1.3 
GeV, m b = 4.3 GeV and m t = 170 GeV (respectively at \i — m c ,m b ,m t ). For the 
parameters V and K, related to the quark and gluon condensates, we have taken 
V = (250 MeV) 3 and K = (230 MeV) 4 and we have assumed that the corresponding 
parameters in the mirror sector scale as V'/V = (A'/A) 3 and K'/K = (A'/A) 4 . In 
figure is shown that A' can be significantly larger than A, allowing for an axion 
mass ~ 1 MeV when f a ~ 10 6 - 10 7 GeV. 

The axion couples with fermions and photons in the standard way, with a 
strength inversely proportional to the PQ constant 0. Referring to ( |1 . 1[ ) we find: 

1 1 X z 1 1 

Cae = "TT Sm 2 /3 , C a d = T7 SIB 2 3 , C au = — COS 2 8 . (2.10) 

For the axion-nucleon couplings we consider c a N ~ 1, while the axion interaction 
with photons is measured by: 

_ 8 6K Tr (M-!Q 2 ) ^ 2z 
ai 3 V + KTr{M~i) " 1 + z' 1 1 

where the trace is taken over the light quark states (u, d), Q are their electric charges 
(+2/3, —1/3), and a = 1/137 is the fine structure constant. In addition, our axion 
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couples with mirror photons, with the constant: 

8 6iT , Tr(M'- 1 g' 2 ) 



on 



3 V + K'Tr(M 



(2.12) 



where the factor C for the case of two (vf, d'), one (u') or no light quarks respectively 
takes the values 2z'/ (1+z'), and 8/3. Hence, the axion decay widths into the visible 
and mirror photons respectively are: 



T(a -> 77) 



64tt 



T(a 



77 J 



9%™l 
64tt 



(2.13) 



where g ai 



' a /2ti fa) and similarly for g' 



In addition, if m a > 2m e , the axion can decay also into an electron-positron pair: 

(2.14) 



V(a 



e e 



8tt \ 



1 - 



Ami 



mi 



where g ae = c ae (m e /f a ). 

We present in figure |3] the axion 
lifetime in the present supersymmet- 
ric model. Notice that, as soon as 
the axion mass turns over 1 MeV its 
lifetime is strongly suppressed. In 
fact in this last case its decay width 
is completely dominated by the a — > 
e + e~ channel. 

For details of the experimental 
bounds we refer to ref. We re- 













tan£ = 25 






tan^ = 55 







10 7 
f a (GeV) 



10 10 



Figure 3: Axion lifetime as a function of f a for 
different values of tan j3 = 5, 25, 55. Supersymme- 
try is supposed to be broken at the scale m s = m t . 



mark here that the parameter region 
of interest for us is f a ~ 10 5 — 10 7 
GeV, a value safe from terrestrial lim- 
its. In this region, the axion mass is 

about 1 MeV so it is not constrained by the standard astrophysical considerations 
applicable for the DFSZ axion. In fact, since it is quite heavy, its production rate in 
the stellar cores, with typical temperatures T up to 10 keV, is suppressed by the ex- 
ponential factor exp(— m a /T). On the other hand, this argument is not applicable to 
the SN, whose core temperature is several 10 keV. If the axion-nucleon couplings are 
large enough, g a p,g an > 10~ 7 (then /„ < 10 7 GeV), the axions are strongly trapped 
in the SN, inside a core of radius R ~ 10 km, and they have a thermal distribution. 
In this case the energy luminosity at t — Is can be estimated as L a ~ fa 6 ^ 11 x 3 • 10 50 
erg/s [p3| , 24| . Hence, for f a < a few 10 6 GeV the axion luminosity L a becomes 



smaller than a few 10 51 erg/s, and the total energy and duration of the SN 1987A 
neutrino burst should not be affected. 
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The constraint related to the axion decay into ordinary photons is not relevant for 
an axion with m a > 2m e , since then, the total axion decay width is almost completely 
driven by a — > e + e~ (see figure |^). On the other hand, another constraint emerges 
due to axion decay into mirror photons. Since the ordinary matter is transparent 
for the mirror photons, the emission of the latter can lead to the unacceptably fast 
cooling of the supernova core. The decay rate for an axion with energy E into mirror 
photons is (m a /E)T', where V = T(a — > 7'7 ; ) is given by eq.( |2.13 ). Therefore, 
we have L'^t) = 47rm a r7 ( f r 2 n a (r,t)dr, where n a = 1.2T 3 /n 2 is the axion number 
density. Taking a core temperature T of about 20 — 30 MeV, one can roughly estimate 
that L' ~ r'm a (47r J R 3 /3)(1.2T 3 / 



7T 2 ) ~ 1Q~ 2 g^m*. Hence, the condition L'^ < 10 51 
erg/s implies that the decay width V, for m a ~ 1 MeV, should not exceed a few s -1 . 
In other terms, using eq. (|2.12|) , we roughly obtain the bound m 2 / f a < 10 -7 MeV. 



Finally, a possible cosmological 
problem must be considered. If the 
axion contribution to the energy den- 
sity at the time of the Big Bang Nu- 
cleosynthesis (BBN) is very large, this 
model can disagree with the strongly 
verified prediction on the primordial 
Helium abundance. This energy con- 
tribution is expressed in terms of ef- 
fective number of extra neutrinos SN U 
in figure f|. 

At the present, deuterium (D) 
and A H e data seem to indicate that a 



0.5 



. 4 



0.3 



0.2 




Figure 4: Effective extra neutrino added due 
to mirror axion as function of axion mass in 
MeV. 



large number of neutrino species is disfavored. There are several recent papers on this 



problem (see, e.g. [25| and reference therein), which give different bounds 5N V < 0.3, 



8N V < 0.5 or even 8N V < 1 [25]. These limits translate to different lower bounds on 
the axion mass. 



3. Mirror Axion: Supernovae and Gamma Ray Bursts 

The Gamma-Ray Bursts (GRBs) puzzle theorists from many points of view [2c]. The 
most striking feature is that an enormous energy, up to io 53-54 erg, is released in a few 
seconds, in terms of photons with typical energies of several hundred keV. The time- 
structure of the prompt emission and the afterglow observations well agree with the 



fireball model |27j in which the GRB originates from the e + e~ plasma that expands 
at ultrarelativistic velocities undergoing internal and external shocks. The Lorentz 
factor of the plasma has to be very large, T ~ 10 2 , which requires a very efficient 
acceleration mechanism. Namely, the e + e~ plasma should not be contaminated by 
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more massive matter (baryons), and hence the fireball has to be formed in a region 
of low baryonic density. 

In particular, the fireball could be powered via annihilation vv — > e + e~ of the 
neutrinos emitted from the dense and hot medium in the accretion disk around a 
central black hole (BH), which can be formed at the merger of two neutron stars 
(NS), or a NS and a black hole (BH) One can consider also the merger of a 

BH and a white dwarf (WD). In addition, the accretion disk can be formed by the 
collapse of a rotating massive star, so called failed supernova or collapsar p9| . These 
objects could potentially provide the necessary energy budget for the GRB. Namely, 
the typical values of the mass M accreted through a disk, the radial size of a disk R 
and the accretion time t, can be estimated as: 



NS + NS 
NS + BH 
Collapsar 
WD + BH 



M~O.1M iZ^SOkm, t ~ 0.1 s 

M~O.5M .R^SOkm, t ~ 0.1 s 

M ~ 2M Q R ~ 200 km, t ~ 20 s 

M~1M Q i?~10 4 km, t~100s. (3.1) 



However, the problem remains how to transform efficiently enough the available 
energy into the powerful GRBs. Due to the low efficiency of the vv — > e + e~ reaction, 
the models invoking it as a source for the GRBs have serious difficulties in reaching 
such large photon luminosities. Even though during the collapse of compact objects 
an energy of ~ 10 53 erg is normally emitted in terms of neutrinos, they deposit only 
a small percent of their energy to fireball, and take the rest away. In addition, 
neutrino annihilation is effective only at small distances, less than 100 km, which are 
still contaminated by baryon load, and cannot provide a sufficiently large Lorentz 



factor 128 



In pi, it was proposed that the emission of the light pseudoscalar particles like 
axions — which can be effectively produced inside the accretion disks and then decay 
into e + e~ outside the system — can provide an extremely efficient mechanism for 
transferring the gravitational energy of the collapsing system into the ultrarelativistic 
e + e~ plasma. The advantage of using the decay a — > e + e~ instead of vv — > e + e~ 
annihilation is obvious. First, it is 100 percent efficient, since the decaying axions 
deposit their energy and momentum entirely into the e + e~ plasma. And second, the 
decay can take place in baryon free zones, at distances of 1000 km or more, and so 
the plasma can get a Lorentz factor T ~ 10 2 . 

Here we suggest that the described mirror axions can be the required pseu- 
doscalar particles. In the dense and hot medium of the accretion disk it is mainly 
produced by its bremsstrahlung in the nucleon-nucleon scattering ||. In order to 
have efficient production, its mass should not over-exceed the characteristic temper- 
ature of the matter, typically a few MeV. Assuming for simplicity the non-degenerate 
and symmetric baryonic matter, the energy-loss rate per unit mass due to emission 
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of axions is: 

e ~ g 2 N pnT^ Y x 2 • 10 31 erg g'V 1 , (3.2) 

where pn is the disk density in units of 10 11 g/cm 3 and T Me v is the matter tempera- 
ture in MeV. 

Axions are in the free streaming regime, if their mean free path is larger than 
the accretion disk size (-R100 = -R/lOOkm), which yields |LJ 



g N < g%(p, T) = 2 x 10" Vi^mCv^oo 72 - (3.3) 

Then the total axion luminosity from the accretion disk with a mass M can be 
roughly estimated as: 

L ~ eM ~ (10 6 g N ) 2 p n T^ v (M/M Q ) x 4 • 10 54 erg s" 1 . (3.4) 

This value can be so big that the axions can extract all the available energy from 
the collapsing system with very high efficiency 

We need the emitted axions to decay into e + e _ outside the disk, in the regions 
of low baryon density, which correspond to distances of several hundred or thousand 
km. 

For the sake of simplicity, let us fix the parameters as m a =1.5 MeV, g ae = 10~ 9 
and g ai ~ 10~ 12 MeV -1 . In this range of parameters the axion lifetime can be 
well approximated by r ~ r(a — > e + e~) = 8Tcg~^(m^ — 4wig)~ 1 / 2 , since the axion 
decay width into electrons is much larger than the one into photons (cfr j2.13|j2.14| ). 



Therefore, since the mean decay length of the axion is D = cTE/m a , where E ~ 2T 
is the average energy of the emitted axions, we find: 

D ~ (10 9 s e r 2 mM 2 eV £MeV x 5 • 10 3 km, (3.5) 

where m MeV represents m a /MeV and E Me y = E/MeV. 

In the view of our mechanism, the short GRBs (duration ~ 10 _1 s), can be 
naturally explained by the NS-NS merger, with typical values M ~ 0.1 M & and 
t ~ 0.1 s, the typical density pu ~ 1 and temperature T ~ 4 MeV. Then the 
total energy emitted in axions from the accretion torus can be roughly estimated as 
£ ~ Lt ~ 10 53 erg, while the mean decay length is D ~ 4 • 10 3 km, much larger than 
the size of the system (R ~ 50 km). Thus, the axions decay in the baryon clean 
zones and deposit their energy entirely to the e + e~ plasma, which can get a large 
Lorentz factor and give rise to rather isotropic photon emission with total energies 
up to 10 53 erg. The relative hardness of the photon spectrum in the short GRBs well 
agrees with this situation. Somewhat more energetic short bursts can be obtained in 
the case of the NS-BH merger, with M ~ O.5M . 
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On the other hand, the long GRBs can be related to the collapsar. In this case, 
we need to estimate the fraction of energy deposited from the accretion disk with 
say M ~ 2M , R ~ 200 km and t ~ 20, within a cone of the solid angle Q along 
the polar axis, where the baryon density is lower and the funneling of the plasma 
in this direction can take place producing a jet expanding outwards. By taking 
p ~ 10 10 g/cm 3 and T ~ 2 MeV, we obtain the beamed GRB with the energy release 
S/Q ~ Lt ~ 10 54 erg. As far as axions decay at large distances, about 1000 km, the 
Lorentz factor can approach large values, V ~ 10 2 . This analysis is supported by the 
result of the simulation in f30fl , which shows that if the energy would be transferred 
to the plasma at distances ~ 600 km, a successful burst could be obtained with 
r ~ 40. 

The axions can be produced also at the supernova explosion. For ~ 10~ 6 they 
are in the trapping regime in the collapsing core and are emitted from the axiosphere 
having a thermal spectrum with a temperature T of a few MeV ||24|| . Therefore, in 
total, an energy of a few xlO 51 erg can be emitted during the collapse period and 
subsequent cooling of the proto-neutron star, in terms of axions with the mean energy 
E ~ 3T. The latter undergo the decay into e + e~ at the distance D ~ 10 3 km. In 
this case, the impact of the axion emission crucially depends on the geometrical size 
of the collapsing star. 

In particular, supernovae type Ib/c result from the core collapse of relatively 
small stars, where the hydrogen and perhaps also the helium shells are missing. 
Their radius can be as small as R ~ 10 4 km, comparable to the axion decay length 
D. This in turn implies that exp(—R/D a ) is not very small, and it can be of order 
10~ 3 to 1, in which case a reasonable amount of axions can decay outside the mantle 
producing a fireball. So the weaker GRBs associated with a supernova type Ib/c 
could take place, having typical energies up to a few 10 51 erg. 

The SN type II are associated with large stars, having an extended hydrogen shell 
(R > 10 7 km). Thus, the axion decay essentially takes place completely inside the 
mantle - the fraction of axions decaying outside the star, exp(—R/D), is essentially 
zero and thus no GRB can be observed. Indeed, the SN 1987A event did not show 
any 7 signal. On the other hand, the energy of a few 10 51 erg released by axion decay 
at distances ~ 1000 km can help to solve the painful problem of mantle ejection (in 
the prompt mechanism, shock usually stalls at a distance of a few hundreds km). 

Concluding, the axion emission from the collapsing systems and their subsequent 
conversion into the relativistic plasma via the decay a —>■ e + e~ outside these systems 
could naturally explain a variety of the GRBs. This mechanism suffers no energy 
deficit and it makes more natural the possibility of the plasma acceleration. In 
particular, the short GRBs, with timescale ~ 0.1 s and total energies up to a few 
x 10 53 erg can originate from the NS-NS or NS-BH mergers, while the collapsar could 
give rise to the longer GRBs, with t ~ 10 — 30 s and S/fl up to a few 10 54 erg. The 
events with t ~ 100 s could also be initiated by the BH-WD merger. In later cases, 
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one can expect more baryon dirty fireball. All this well agrees with the observational 
features of the GRBs. As explained, an interesting possibility is the association of 
some weak GRBs (with total energy up to 10 51 erg or so) with the supernovae type 
Ib,c. Considering that the axion emission could also help the supernovae type II 
explosion, we see that this model can provide a unified theoretical base for the GRB 
and SN phenomena. Interestingly, the emission of these axions can also be important 
for explaining the observed GRB's preceded by supernova explosions, via collapse of 
usual neutron stars to quark or hybrid stars [|3l| . 

4. Conclusions 

We have presented a new model of axion for the solution of the strong CP problem. 
The main feature of this model is the modification of the relation between the axion 
mass and the PQ constant, with quite interesting phenomenological consequences. 

We have hypothesized the existence of another sector of particles and inter- 
actions, the mirror world, which is an exact copy of our world but with a larger 
electroweak scale v' ^> v. This difference also implies a different dynamics in the 
QCD sector. The fermion masses are in fact driven by the Higgs VEVs, and the 
different thresholds for the mirror masses lead to a pole A' (in the evolution of the 
mirror strong coupling constant), that can be significantly higher than the ordinary 
one A. Since the value of A in the two sectors depends on the RG evolution equa- 
tions, and consequently on the matter content of the theory, this behavior can be 
improved in a supersymmetric model. 

The result is that the relation of axion mass to the PQ constant is relaxed, 
and our particle can be quite heavy, maintaining the weak coupling with matter 
and photons typical of the invisible axions. This behavior has great benefit on a 
phenomenological and astrophysical ground. Because of its large mass our axion has 
in fact no influence on stellar evolution. All astrophysical constraints come from SN 
explosion and allow its mass to be m a ^ MeV with a PQ constant of order 10 5 — 10 6 
GeV. 

Such a heavy axion happens to have quite an interesting phenomenology, in 
particular in relation to the GRB and SN physics. The axion emission from the col- 
lapsing systems, and their subsequent conversion into the relativistic plasma via the 
decay a — > e + e~ outside these systems, could naturally explain the fireball formation 
and, consequently, a variety of the GRBs. This mechanism was first proposed using 
the reaction vv —* e + e~ , instead of the axion decay a — > e + e~, but the advantage of 
using the last decay, instead of the neutrino annihilation, is clear. First, it is 100 per- 
cent efficient, since the decaying axions deposit their energy and momentum entirely 
into the e + e~ plasma. And second, the decay can take place in baryon-free zones, at 
distances of 1000 km or more, so the plasma can get a Lorentz factor T ~ 10 2 . 
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Another interesting possibility is the association of some weak GRBs (with total 
energy up to 10 51 erg or so) with the supernovae type Ib,c. The axion mean-free 
path is a few 10 3 Km, so a few of them can reach the supernovae type I surface and 
then decay into photons, giving rise to the observed weak GRBs. Finally, the axion 
emission could also help the supernovae type II explosion, thus providing a unified 
theoretical base for the GRB and SN phenomena. 

As a final note, the parameter window allowed for our axion is accessible to the 
axion search in the future reactor and beam dump experiments. 
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